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Abstract 

We consider the finite density, zero-temperature behaviour of quark matter in the 
instanton picture. Since the instanton-induced interactions are attractive in both 
qq and qq channels, a competition ensues between phases of matter with conden- 
sation in either or both. It results in chiral symmetry restoration due to the onset 
of diquark condensation, a 'colour supercondutor', at finite density. Also present 
is a state with both manners of condensation, however this phase is found to be 
thermodynamically disfavoured for equilibrium matter. Properties of quark matter 
in each phase are discussed, with emphasis on the microscopic effects of the effective 
mass and superconducting energy gap. 



PACS: 11.30.Rd; 11.15.Tk; 21.65.+f 



1 Introduction 



The idea that the QCD partition function is dominated by instanton fluctuations of 
the gluon field, with quantum oscillations about them, has successfully confronted 
the majority of facts we know about the zero-temperature, zero-density hadronic 
world (for reviews see refs.fl], [J). Instantons have been reliably identified in lattice 
simulations (for a review see ref.0), and their relevance to hadronic observables 
clearly demonstrated From the theory side, the instanton vacuum constructed 
from the Feynman variational principle [[5], || , gives an example of how the necessary 
'transmutation of dimensions' can actually happen in QCD, meaning that all dimen- 
sional quantities can be expressed through the QCD scale parameter, Aq C d- It can 
be added that in the solvable N = 2 supersymmetric version of QCD it is instantons 
- and nothing besides them — that are sufficient to reproduce the expansion of 
the exact Seiberg-Witten prepotential [0]. 

Of even more relevance to the topic of this paper, the instanton vacuum provides 
a theoretically beautiful and phenomenologically realistic mechanism of spontaneous 
chiral symmetry breaking in QCD (for a review see Jl|). In addition, it has recently 
gained a strong support from direct lattice studies M. Therefore an expansion of 
instanton ideas to a new frontier, in this case to nonzero matter density, seems to 
be well justified. 

The key point is that instantons induce interactions that are attractive not only 
in the qq channel (leading to the spontaneous chiral symmetry breaking) but also 
in the qq channel (potentially leading to diquark condensation). This has been first 
noticed in ref. ||. It becomes especially clear in the special case of two colours 
(N c = 2) and two light quark flavours (Nf = 2). In this case instanton-induced 
interactions possess not a global SU (2) x SU (2) as would be for any N c > 2, but 
a larger S77(4) symmetry [Kj (often referred to as Pauli-Giirsey symmetry). At 
N c = 2 the attraction in the qq and qq channels are exactly equal. For this reason 
qq and qq condensates belong in fact to one phase: one condensate can be rotated 
to another along the Goldstone valley [p~Qj|p|. 



Switching in a nonzero chemical potential /i violates explicitly the global SU(A) 
symmetry of the N c — 2 world, the degeneracy of the 5-dimensional Goldstone 
valey is lifted, and one can ask which of the condensates becomes preferred. In ref. 
[p~2]j arguments have been given that it is the diquark condensate. In a previous 



1 In the original paper |T(J a general case of the 5'[/(4) breaking has been considered, leading 
to nine Goldstone particles. A closer inspection shows, however, that the instanton-induced inter- 
actions possess an additional degeneracy, leading to only five Goldstone particles, instead of nine 
which would be the general case. Instantons apparently 'know' about the Vafa-Witten theorem 
which guarantees that the symmetry breaking pattern at N c — 2 corresponds to five Goldstone 
particles O]. 
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publication []T3| we have confirmed this expectation by direct calculations; see also 
below. It should be noted that at N c = 2 singlet diquarks are nothing but the 
colourless 'baryons' which happen to be bosons, and their condensation does not 
break the colour symmetry. 

In the opposite limit, N c — > oo, the diquark interaction is suppressed by a factor 
of ~ 1/N C as compared to the qq one, and so is the diquark condensate. The question 
is then, is our world with N c = 3 closer to N c = 2 or iV c = oo? 

The possibility of diquark condensation at any N c , as induced by instantons, has 
been studied in ref. [Q. Contrary to the case of N c = 2, a diquark condensation 
at N c > 3 would inevitably spontaneously break the colour symmetry, similarly to 
the Higgs breaking of the SU(2) gauge symmetry of electroweak interactions. For 
that reason the (possible) symmetry breaking of the colour SU(N C ) has been named 



'dynamical Higgs mechanism' in ref.[|14[], and a parallel with the superconductivity 
has been drawn. The presently used term is 'colour superconductivity' to which we 
shall adhere. 

For zero chemical potential only a metastable diquark-condensed state has been 
found at N c > 3 |TJ|]. Furthermore, at N c = 3 the scalar diquark already appears 
to be unbound in the vacuum, indicating that our world is in a sense closer to the 
idealized N c — » oo than to the N c = 2 world. Parametrically, the diquark mass 
is ~ 1/p ~ 1 GeV where p is the average instanton size as explained below. This 
means that a scalar diquark correlation function should decay with the exponent 
corresponding to the 'constituent' quark threshold 2M(0) ~ 700 MeV, which seems 
to be supported by recent lattice measurements [|l^]. Nevertheless, it has been 



suggested in Jl4| that qq condensates could be found as metastable states in heavy 
ion collisions and in astrophysics. 

More recently it has been argued by the Princeton/MIT and Stony Brook groups 



1^, [16|] both using instantons as a framework, that taking nonzero fermion density 
shifts the balance in favour of the diquark condensation, and that at certain critical 
chemical potential fi c there should be a phase transition from the usual broken 
chiral phase to the colour superconducting stateQ Avoiding some of the unnecessary 



approximations made in Jl2], [16| we arrive essentially at the same conclusions ]I3] . 

In this paper we study the competition between various quark channels in a 
more systematic way than in our previous publication |Tj|. Since 'colour supercon- 



ductivity' implies colour symmetry is broken, one can imagine several phases with 
chiral symmetry broken or restored for quarks of different colours. We explore these 
possibilities using what amounts to a virial expansion in the instanton density, and 
carrying out detailed calculations to first order. 



2 See the Proceedings of the International Workshop of QCD at Finite Baryon Density (Biele- 
feld, Germany), where these matters have been extensively discussed [ fl7|| . 
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The chiral broken phase is characterized by a nonzero order parameter (qq) ^ 
or, equivalently, by a nonzero dynamical or 'constituent' quark mass which in fact 
is momentum- dependent. We find that, as far as one remains inside this phase, the 
quark occupation number n(p) is a perfect Fermi step function, despite strong energy 
and momentum-dependent interactions between quarks. It should be contrasted to 
the colour superconducting phase with (qq) ^ but (qq) = 0. In this case the quark 
occupation number n(p) is distorted near the Fermi surface and is of a typical BCS 
type. 

In our numerics we use the standard characteristics of the instanton ensemble. 
They lead to very reasonable values of the quark condensate and of the constituent 
quark mass at zero chemical potential. However, the same instanton characteristics 
lead inevitably to a very early transition to the superconducting state: it occurs at 
quark densities less than the normal nuclear density; after the system jumps into the 
superconducting phase the density appears to be about twice the nuclear density. 
Taken literally, it suggests that the ordinary nuclear matter is a 'boiling' mixture of 
the two phases , however, it most probably contradicts the lore of conventional 
nuclear theory^ To 'save' the nuclear matter, one would probably need to go beyond 
the mean-field approximation actually used in this paper. 



2 QCD Instanton Vacuum 

The construction the QCD instanton vacuum has been reviewed in previous publi- 
cations 0, 13], therefore here we simply review the main steps. 

• The use of the Feynman variational principle applied to instantons leads to the 
stabilization of the grand canonical instanton ensemble, with the main charac- 
teristics of the instanton medium, namely, the average 4-dimensional instanton 
density, N/V = 1/i? 4 , and the average instanton size, p, expressed through 
the only dimensional quantity Aqcd- The 2-loop calculations performed in 
|, | give 

_ /= 0.48 - 1.35 

IV MS IV MS 

Taking = 280 MeV one finds p « 0.35 fm, R « 0.95 fm, p/R » 

1/3. This small ratio has been previously suggested on phenomenological 
grounds by Shuryak fl9 |. The smallness of the p/R ratio implies that the 



3 We take an opportunity to thank G.Brown and K.Langanke for a instructive discussion of this 
point. 
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packing fraction of instantons in the vacuum, i.e. the fraction of the 4d volume 
occupied by the balls of radius p, is quite small: 

1 2 9 V 2 R 4 ~ 10' U 

Though the packing fraction (||D is but numerically small, one can treat it as 
a formal algebraic parameter, and develop a perturbation theory in it. Taking 
nonzero matter density and/or nonzero temperature can only decrease the 
ratio p/R. 

When one switches in light quarks on top of the instanton ensemble, the 
existence of the small packing fraction parameter enables one to separate con- 
tributions of high and low fermion eigenmodes. The high-frequency part is 
controllably factorizable into contributions from individual instantons. This 
part can be seen to renormalize somewhat the one-instanton weight or in- 
stanton 'fugacity'. The low-frequency part is dominated by the would-be zero 
modes of individual instantons. The normalized (that is localized) zero modes 
in the background of individual instantons exist both in the vacuum (p = 0) 
p0| and at p ^ |3T], |22| , see the Appendix. In the ensemble of instantons and 
anti-instantons they cease to be exactly zero modes (that is why we call them 
the 'would-be zero modes'). The spectral density of the Dirac operator at 
small eigenvalues is obtained through the diagonalization of the matrix made 
of the overlaps of the would-be zero modes belonging to individual instantons 



and anti-instantons p3fl . A nonzero spectral density at zero Dirac eigenvalue 
signals chiral symmetry breaking: it is actually due to the derealization of 
the would-be zero modes owing to the 'hopping' of quarks from one instanton 
to another E3 . 



In fact, three seemingly different but actually equivalent methods have been 
developed to describe chiral symmetry breaking by instantons at zero p: i) 
diagonalization of the random matrix made of the would-be zero modes' over- 
laps p3 , |25| ; ii) finding the quark propagator in the random instanton ensem- 



ble p4| , |26j ) which exhibits the appearance of a momentum-dependent quark 
mass; Hi) derivation of the effective low- momentum theory for quarks with 
instanton-induced interactions [[[], p7[| . 

The three approaches underline different sides of the physics involved though 
mathematically they prove to be equivalent. It is straightforward to generalize 
them to the case of nonzero chemical potential. In this paper we use only the third 
method, it being the most economical. 
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3 Effective Fermion Action for Finite Density 



With the instanton solution representing a dominant, stable fluctuation of the gauge 
field, the original strategy for computing quark observables in the QCD vacuum was 
to calculate the quantity in the presence of a random instanton background and then 
average over an ensemble of such field configurations. However, the analysis here 
calls for a more transparent manifestation of the effects from instantons in order 
to resolve the possible mechanisms for symmetry breaking. This is obtained by 
first averaging over the ensemble of instanton and anti-instanton configurations to 
formulate an effective theory in terms of interacting quarks, where explicit instanton 
effects are absorbed into the form of the interaction. 

This evidence of averaged instantons is retained in the interaction through the 
would-be fermion zero modes, one for each flavour. Thus the consequent interaction 
is a vertex involving 2Nf quarks, commonly cited as a 't Hooft interaction after 



the first author to specify their proper quantum numbers f20|| . For two flavours, 
it may be cast to resemble the Vaks-Larkin/Nambu-Jona-Lasinio model which has 
long been a popular model for quark and hadron phenomenology. 

The field of a single instanton/anti- instanton determines a particular solution 
for the gauge field, Aj, and in the absence of quantum fluctuations the zero mode, 
is simply the solution of the Dirac equation with zero eigenvalue: 

(i@ - 1^74 - 4j) ft) = • fi) = 0. (3) 



The exact solution for nonzero /i was found by Abrikosov [|2T| and is written in the 



Appendix. It can be decomposed into chiral compontents as 



(4) 



meaning that a (right) left handed zero mode is generated by an (anti-) instanton. 
This is a result of the instanton solution's structure and the self-dual equations from 
which it comes. With $jm we will denote the conjugate zero mode, which takes 
the chemical potential argument of opposite sign: = — //), where 

the dagger means hermitean conjugate. This is a direct consequence of the non- 
hermiticity of the /i-dependence which arises in the Dirac operator. 

This solution is not physically realizable, since the QCD background is not mod- 
elled by one but many (a 'liquid' of) instantons. Localized around each there is but 
a would-be zero mode, from which the effective action is constructed. This action 
is determined by building a low-momentum partition function, the foundation of 
which is the Green function for a quark in the field of one instanton. It can be 
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expressed as a sum over the complete set of eigenfunctions, and approximated as 
S (x,y) = ^ if 3/ = — —■ 

+ S {x,y;ti) . (5) 



im 

At low momenta this propagator is dominated by the zero modes while at large 
momenta it is essentially reduced to the free one, So(x,y; p). Therefore, eq. fl5|) can 
be regarded as an interpolation of the true propagator in the field of one instanton. 
It becomes exact both at large, p ^> 1/p, and small, p -C 1/p, quark momenta. 

The partition function which produces the necessary propagator (^) is of the 
form 

N+ „_ \ N- 

(6) 



V \ V 



Z = j Dip Dip^ exp J ipf(i$ — ifj,j^)ipf 

The pre-exponential factors contain the instanton-induced interactions between fermions, 
and for Nf flavours are the specific nonlocal 2iVj-fermion vertices 

Nf 

Y£ f [ij>, ^} = (-) N fjd 4 z du]j J d ' x d 'y 

•1>lA x )( id - ^) T ^iA x - z )®i,i(y - z )( id - w)^l,r{v)- ( 7 ) 

Here we use the notation x^ = x^a^, where the 2x2 matrices cr^ = (±i<7, 1) 
decompose the Dirac matrices into chiral components, and it is understood that p 
written as a four- vector is p a = (0, p). Note that the zero mode in the field of an 
(anti) instanton couples to that of a (right-) left-handed quark, and the effective 
range of the interaction is the average instanton size, p. 

It should be stressed that correlations between instantons induced by fermions 
are inherent in this approach; as to correlations induced by gluons, they are effec- 
tively taken care of by the use of the variational principle [f| |(| resulting in the 
effective size distribution. For simplicity we freeze all the sizes at the average value 
p, but average explicitly over random position (z) and orientation (U) variables. 

Fermion operators in the pre-exponent are not convenient; these operators can 
be raised into the exponent with the help of a supplementary integration over a pair 
of Lagrange multipliers, denoted A±: 

Z = J d\ + d\_ J Dip Dfl exp j J d 4 x - 2/174)^ + A+Y^ + X-Y Nf 
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Indeed, integrating over A± by the saddle-point method one recovers eq. (|]). No- 
tice that the saddle-point integration becomes exact in the thermodynamic limit: 
N±, V — > oo with N/V fixed. Through this procedure we obtain a purely exponen- 
tial integrand which is the required effective action. 

Two important consequences follow and should be emphasized. First, the cou- 
pling constant of the 2A^-fermion interaction, whose role is played by A±, is not fixed 
once and forever; its value is found from minimizing the free energy after integration 
over fermions is performed. Therefore, the strength of the interaction depends itself 
on the phase the fermion system assumes. Second, the saddle-point value of A± is 
not proportional to the instanton density, as one might naively think. For exam- 
ple, in the case of chiral symmetry breaking it behaves as A ~ (N/V) 1 ^ Nf ^ 2 [p7 



Both circumstances are due to the pecularity of the instanton-induced interactions 
dominated by the existence of the zero modes. Would these modes remain zero, 
the fermion determinant would be zero, thus suppressing the presence of instantons 
themselves. It is this intricate fermion self-supporting mechanism which makes the 
instanton-induced interactions different from a more naive NJL model where the 
strength of the effective four-fermion interactions are chosen once and forever. As 
will be shown below, this is the mechanism by which chiral symmetry is restored at 
large chemical potential in the case of two flavours. 

For practical applications it is favourable to use the Fourier-transformed ex- 
pressions of the quark zero modes, which are written explicitly in the Appendix. 
These complex functions of the four-momenta and chemical potential determine the 
(matrix) formfactors attached to each fermion leg of the vertex: 

F(p,fj) = (P+ iii)~ip(p,n) + , J^{p, -fi) =<p*(p, -^T(p + i/i) + . (9) 
With such definitions, the interaction terms may be written in momentum space, 

Xf r (d 4 Pf d% ' 



IV / 

F + [^ f ] = fdu /n 



(2*0* 



(E(p/-*y 



•n[v>u,v(p/)^(p/>^ ( io ) 

with a similar form for Y~ which carries right-handed quarks. The first indices on 
the fermion operators refer to flavour, the Greek to colour (1 . . .N c ), and the last 
denote spin (1,2). This formulation of the effective interaction retains the full p 
and /i dependence of the zero modes, as opposed to the approximate treatments in 



other recent works |T2|, [L6|. At the same time this formalism enables one to reduce 
the problem of determining the phase structure to algebraic equations where the 
coefficients are given by certain integrals over the formfactor functions (0). 
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4 One Light Flavour (A Pedagogical Aside) 



In the case of one flavour, the effective action assumes a relatively simple and hence 
instructive form. The formfactors combine and appear as an effective mass term. 
The interaction remains non-local, thus one finds an induced mass, A4, with a 
specific momentum dependence: 

Sint = -A+ J "^j J 'dU ^{jp)T{p,^ k U^e ml ^{p, -a*)™V?(p) 
-A_(L <-> R) 

= M(p,n)^(pW(p), (11) 

where 

M(p,n) = X(p + ifi) a (j> + ifi) a (pfi(p,fi)(pi3(p,fi). (12) 

At finite // the formfactors are complex, since the zero modes are solutions of eq. (13). 
This interaction is sufficient to spontaenously break chiral symmetry. However, the 
coupling constant A has here been introduced as a Lagrange multiplier in order to 
obtain a manageable form for the effective action, and is in principle an integration 
over all possible coupling strengths. This simplifies to a saddle-point approximation 
which becomes exact in the thermodynamic limit of N, V — > oo. Such an evalution 
of the integral also naturally connects the instanton density to this coupling, in that 
the resulting 'gap' equation is: 

— = 47V / M ^) 2 (13) 

This equation is a direct generalization of the gap equation found for ji = in ref. 
24jj . Although Ai and the [i term are complex, one finds that the imaginary part 



of the integrand is odd in p 4 and hence vanishes under energy integration. This is 
equivalent to the statement that while individual eigenvalues of the Dirac operator 
are imaginary, their product remains real. By preserving the exact form of the zero 
mode solution, this property becomes manifest in our treatment. Through solving 
this equation self-consistently, one obtains a solution for the A at any given quark 
chemical potential fi. This in turn is proportional to an effective mass evaluated in 
the zero momentum limit, |.M(0,/i)| = (2iip) 2 X/N c . 

To first order in A, this procedure amounts to averaging the effects of every 
instanton and anti-instanton which modify the quark propagator. In principle, the 
instanton density N/V is modified by a finite quark density. It scales with the 
fugacity, x = as X 4 ^ \ where b is the Gell-Mann-Low coefficient of 

QCD. Fugacity not equal to unity simply arises from the change with /i of the 
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fermion determinant in the partition function. Upon a systematic expansion of 
this determinant in A one finds the first order term vanishes, for this is exactly the 
saddle-point condition of eq. flTB|). Thus modifications of the instanton background 
appear only when the determinant is expanded to order A 2 , which is beyond the 
scope of this work. 

Without a first-order effect on N/V, there exists no mechanism in the theory to 
surpress the instanton background and modify the form of the gap equation above. 
While a finite quark chemical potential can and does modify the coupling strength 
and hence the effective mass, the coupling is never forced to zero. Therefore chiral 
symmetry remains broken at any finite baryon density. This result can of course 
change when effects of higher order in the instanton density are taken into account. 
The picture is also radically altered with more than one flavour, as possible quark 
pairing offers an alternative mechanism for chiral symmetry restoration. 



5 Two Light Flavours 

The remainder of this paper investigates the case of massless Nf = 2, corresponding 
to a system of chiral up and down quarks. Furthermore, we will assume the CP 
invariant case of 9 = 0, which requires N + = AL = N/2 and hence A+ = A_ = A. 
With the definitions and notation of eq. (|9]) and eq. fllUP , the interaction terms may 
be written in momentum space, 

\Y + = XJ (2tt) 6 {p 1 +p 2 -h- k 2 ) 

with a similar form for XY~ . 



5.1 Gap Equations 

Since the instanton-induced interactions ( |14D support both qq and qq condensation, 



it is necessary to consider the two competing channels simultaneously. This means 
that one must calculate both the normal (S) and anomalous (F) quark Green func- 
tions. A colour /flavour /spin ansatz compatible with the possibility of chiral and 
colour symmetry breaking is 

<^(p)^ ft .(p)> = SffiSfa)) for a,P = l,2, 
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(ip fai (p)4pM = H S P S M for a, (3 > 2 , 
(^(P)^f(-P)) = (4 ffl (p)< J (-p)) = ef'fPM&Ftp) , (15) 

where [7] refers to some generalized direction(s) in colour space, and it is this set 
of N c — 2 indices which signals the breaking of colour symmetry. In the particular 
case of N c = 3, where the colour symmetry is broken as SU (3) — > SU (2) x £7(1) and 
our ansatz considers the 3 channel, we will by convention take [7] = 3; for N c = 4 
one can take [7] = 34 and so forth. In the event of colour symmetry breaking, the 
standard propagators (and ensuing condensates) will lose their colour degeneracy 
and the separation of S(p) into S\(p) and S 2 (p) becomes necessary; otherwise the 
Schwinger-Dyson-Gorkov equations do not close. 

Written in the chiral L, R basis, the 4x4 propagators Sx j2 (p) are °f the form: 



S(p) 



G(p)i z(p)s ( P y 

Z(p)S (p)- G(p)l 



(16) 



Here the off-diagonal, bare propagator So(p) ± = [(p + z/i) ± ] _1 is modified by the 
scalar functions Zi j2 (p), and is augmented on the diagonal by the scalar G\^{p) 
which if nonzero break chiral symmetry. 

Using the instanton-induced interaction (|14"D one can build a systematic expan- 
sion for the F,G Green functions in the 1/N C and p/R parameters. In the leading 
order in both parameters we restrict ourselves to the one-loop approximation shown 
in Fig. |I| It corresponds to a set of self-consistent Schwinger-Dyson-Gorkov equa- 
tions. An important //-dependence enters through the formfactors in flUp . With 
F(p) = F*(p) these diagrams lead to the set of five algebraic equations for the scalar 
Green functions Z lj2 , Gi,2 and F: 

Z x (p) = 1 - G 1 (p)A(p, f i)M 1 - 2F(p)B(p, fi) A 
Z 2 (p) = l-G 2 {p)A{p,v)M 2 
Gi(p) = Z x {p)ip a {p, n)tp a (p, n)M x 
G 2 (p) = Z 2 (p)ip a (p, fi)ip a (p, n)M 2 

F(p) = 2Z 1 (-p)tp a {p,fi)ip a (-p,ij,)A. (17) 
The constants Mi, M 2 , and A will be defined below and the functions 

B(p, y) = (p 2 + V 2 )<P/](p, fj.)ipp(-p, fi) + (p + ifi) a (p a (p, p) (p - in) p ^p{-p, p) 

-(p + ifi) a if a (-p, fi)(p - ifi)0tpp(p, fi) , (18) 
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Figure 1: Schwinger- Dyson- Gorkov diagrams to first order in A, corresponding to five 
scalar equations in the text. The first is for the off-diagonal and the second for the 
diagonal components of Si and S 2 , the third is the anomalous diquark propagator. 

are the formfactors which arise from the zero modes (see the Appendix). At fx = 
we have A(p, 0) = B(p, 0), but for any finite \x the direction of the momentum flow 
through each vertex leg is critical. 

The condensates g\, g 2 , and / are the closed loops contributing to the quark self- 
energy. They are found by integrating the appropriate Green function, modified by 
the vertex formfactors, over an independent momentum: 

! = Nhl§r B ^ F ^- (19 > 

The quantities Mi j2 and A are linear combinations of the condensates g± i2 and /: 

m ' = ( 5 -£k + ( 2JV <- 5+ J:K 

M 2 = 2 (2 - -i) 9 i + 2(JV„ - 2) m , 

A =(! + £)/. (20) 

The Mi 5 2 are measures of chiral symmetry breaking, which act as an effective mass 
modifying the standard quark propagation. On the other hand the diquark loop 2A 
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plays the role of twice the single-quark energy gap formed around the Fermi surface. 
The Fermi momentum, in the absence of chiral symmetry breaking, will remain at 
Pf = ji regardless of the magnitude of A. It should be noted that the effective 
masses M 12 exhibit different behaviour in N c than the superconducting gap A: the 
former are N c times larger than the latter. Although one can of course express a 
solution in terms of (#1, <72> /) as one in (M 1; M 2 , A), we will retain both notations 
in the following text. The reason is the distinct physical interpretations for the two 
variable sets: the first quantifies the symmetry breaking by channel while the second 
measures its effects on the quark states with particular colour. 

After determining the five scalar functions through solving eqs. ( |P?D and insert- 
ing these solutions into eqs. ( |19D we find the coupled equations for the condensates 
themselves: 

A Mi f d 4 k a(p, fi) [1 - 4/3(p, n)A 2 + a*(p, n)M?} 

9i 

92 



Nl-lJ (2tt) 4 [1 + a{p,n)Ml] [1 + a*(p,n)M?} - 16/3(p, /i) 2 A 4 
AM 2 f d A k a(p, fi) 



Nl-lJ (2tt) 4 l + a(p,n)M$ ' 
f 2AA r d'k gfa /i) [1 - 4/3(p, /x)A 2 + q(p, 

1 N^-lJ (2vr) 4 [1 + a(p, n)M%\ [1 + a*(p, fi)M?} - 16/3(p, /i) 2 A 4 ' 1 ] 
where yet another pair of functions has been introduced, 

a(p,fi) = A(p, fi)<p a (p, n)<p a (p, fi) , Pip, ft) = B(p,n)<p a (p,n)<p a (-p,fi) . (22) 

Note that while (3 is real, the function a is complex. Although the integrands above 
are complex, all imaginary parts are odd in p^ and thus vanish under integration. 
As usual for gap equations, there is the possibility of a solution where some (or even 
all) of the condensates vanish. 

The magnitudes of these condensates are not yet determined, as this requires fix- 
ing the coupling constant A. This is done through minimizing the partition function 
(H) and leads to the saddle point condition on A, replacing the integration with 

y = \(Y+ + Y-}. (23) 

On the left-hand side is the instanton density, which we will here take to be fixed at 
its vacuum value, although in principle it will have some correction due to the finite 
quark density. Evaluating the right-hand side requires calculating the one-vertex 
contributions to the free energy, which in this case includes two-loop, figure-eight 
type diagrams formed by joining the four fermion legs into two pairs as shown in 
Fig. |2|. This can be written concisely in terms of the condensates: 

- = A(F+ + Y-) = -A-^ 1 \2g x M x + (N c - 2)g 2 M 2 + Af A] . (24) 
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n=(S X S) + (F 




Figure 2: Contributing diagrams to Q at order A 1 . 

As with the single flavour case we are quenching by truncating our perturbative 
treatment in A, and thus the instanton denisty on the left is held fixed. This equation 
for A and the definitions of the condensates comprise a system of equations 



which can be solved for all quantities. It follows from eqs. (f20|) , (|21 ) and (gj) 
that the mass and energy gaps and A are proportional to the square root 
of the instanton density, N/V, while the saddle-point value of A is, to the first 
approximation, independent of N/V. 

Once these quantities are found the chiral condensate may be calculated, being 
a trace over the quark propagator. This is distinct from the effective masses M 12 , 
as it does not contain contributions from the formfactors present in the loops from 
which the effective masses are computed. Present instead is the solution for the 
propagator, in particular its diagonal elements: 

- (^1p)Mink = i{^H)Eucl = i 1 7^J Tr [Sip)] 



(2tt)<- 
d 4 p 



[2G 1 (p) + (N c -2)G 2 (p)). (25) 



As with the previous integrals, the integrand reflects in p<± such that the result is 
guaranteed real. 



5.2 Thermodynamic Competition 

As a consequence of having two different modes of quark condensation, one obtains 
multiple solutions for the characteristics of the quark medium at any fixed chemical 
potential. Specifically, there will be a competition between the following phases: 

(0) Free massless quarks: g\ — g 2 = f = 0. 

(1) Pure chiral symmetry breaking: g x = g 2 0, / = 0. This is the standard 
vacuum scenario. 

(2) Pure diquark condensation: g\ — g 2 — 0, / ^ 0. This is the 'colour supercon- 
ducting' phase, and leaves quarks of colour [7] free, since they neither participate 
in the diquark formation nor exhibit broken chiral symmetry. 
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(3) Mixed symmetry breaking: g\ ^ g 2 7^ 0, / 7^ 0. In fact one need not have 
finite g\ , since it corresponds to chiral symmetry breaking by the quarks of transverse 
colour to the diquark condensate, however such a phase requires that there at least 
be chiral breaking in the parallel-coloured quarks. 

Phase (0), with all symmetries restored, would mean that the average four- 
fermion vertex, (Y + + Y~), is zero. It would imply, via eq. (|J), that the saddle- 



point value of A is infinite. Being substituted into eq. (|23|) , infinite A means that the 
fermion determinant vanishes. Consequently, the whole instanton vacuum is severely 
supressed, which means a large loss in the free energy. One can put it in another 
way: if the density of instantons N/V is considered fixed, there is no solution to the 



combined eqs. fl20| , |21| , |23[) with all condensates being zero, at least in the one-loop 
approximation to the Schwinger- Dyson- Gorkov equation we are considering. 

To resolve between the remaining, symmetry-breaking solutions the free energy 
is minimized. Consistent with the evaluation of the Green functions, it is calculated 
to first order in A. Repeating the calculation of figure-eight diagrams and recalling 
the explicit dependence on A in eq. (H) we obtain the free energy 



InZ 



V 3 pv 3 

Here V3 is the three-volume while V is the Euclidean four-volume, and Qq is the free 
energy for a gas of free quarks. The last two terms are precisely the quantity which 
must vanish under the saddle-point determination of A, and thus we have 

n=3> + ^ to f*V (27) 

V 3 V, V \N/V) 

Thus the phase which features the lowest coupling A is the thermodynamically 
favoured. 

If we first concentrate on a competition between the two phases of simple chiral or 
colour symmetry breaking, (1) and (2), the former is marked by nonzero g = gi — g 2 
and the latter by nonzero /. Each case leads to its own value of A through solving 
eq. p^ ) for any given /i. For the first, the gap equation reduces to 

\ _ ^ - 1) y 

x -^jjv — ' (28) 



while for the second one finds 



_ 16(iV e 2 -l)(iV c + l)/ 2 
A " NjNjV) ' (29) 



15 



At certain value of p = p c these solutions for A cross; this is the point where the 
phase transition occurs. This point is defined by the condition 



f(p c ) _ I N C (N C - 1) 

SS"V 2 • (30) 

When the ratio between the condensates is less than the constant on the right, 
phase (2) is favoured; otherwise it is chiral symmetry that is spontaneously broken 
in phase (1). However, the phase structure is not so simple in general, considering 
the possible presence of phase (3). Calculations were carried out for N c = 2 and 
N c = 3, taking the values for N/V and p specified in Section 2. From eq. (|3C|) 
one observes that at large N c it becomes increasingly difficult to get to the colour 
superconducting phase. 



6 Two Colours 

In this case it is obvious that colour symmetry is not broken by diquark formation, 
which here correspond to colour-singlet 'baryons', and hence there is only one pos- 
sible chiral condensate g\ — g. Not so obvious is that at p — the colour and 
flavour SU(2) groups are arranged into the higher S77(4) symmetry. So no mixed 
phases exist; the theory has only one symmetry. The instanton vacuum accounts 
for this [10 1, and in the context of the analysis here this corresponds to f 2 + g 2 
being the only discernable quantity in the gap equations (pip. Numerically, we find 
VP + g 2 = 147 MeV. 

At finite p, however, the 577(4) symmetry is explicity broken and there is a 
thermodynamic competition between phases of pure chiral or diquark condensation. 
For all finite chemical potential, numerical calculations reveal that / decreases faster 
than g. Since in this case the critical ratio of eq. (|30| ) is unity, we conclude that for 
any finite density the N c = 2 world prefers diquark condensation to chiral symmetry 



breaking. This finding is in agreement with the reasoning of ref. [O] and the lattice 



results of ref. 28 



7 Three Colours 

With three colours, the more relevant case for QCD, the phase structure is richer 
due to the possible appearance of phase (3). Furthermore, the physics is naturally 
completely different, since here any diquark formation will spontaneously prefer a 
particular colour. This dynamical consequence mimics the Higgs mechanism as a 
spontaneous breaking of the SU(N C ) gauge group. 
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Figure 3: Chiral and diquark gaps as a function of chemical potential in the phase 
of mixed symmetry breaking, (3). The dashed line is the phase (2) solution of the 
diquark energy gap. 

7.1 Numerical Results 

The system of gap equations is comprised of all three in eq. (EH]), which are ac- 
companied by the condition (|2"4j). For a given /i a set of solutions can be obtained 
numerically, each of which corresponds to a phase of the quark matter characterized 
by the mass or energy gaps formed. 

In the vacuum, where fi strictly vanishes, there exist solutions for all three 
symmetry-breaking phases. The thermodynamically favoured is one of spontaneous 
chiral symmetry breaking, and hence the equilibrium description of the vacuum 
in the instanton picture is recovered. Numerical calculations show that the con- 
stituent quark mass is M = Mi = M 2 = 346 MeV, corresponding to A (1) = 0.311. 
Two metastable states are also present, these being manifestations of phases (2) 
and (3). They feature = 0.674 and A^ 3 ) = 0.673, and are thus clearly only local 
minima of the free energy. It is noteworthy that the diquark gap in the metastable 
state is quite large: A = 220 MeV. This is in agreement with that found for fi = 
inref.0g 

Although it is an unstable state in vacuum and remains so as fi is increased, the 
phase of mixed symmetry breaking is worthy of attention. It is marked by nonzero 

4 In that work the notation m c — 2A was used. 
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condensates g 2 and /, which leads to not only a superconductivity energy gap but 
also finite and distinct effective quark masses, M\ and M 2 . The magnitudes of 
these gaps are plotted as a function of chemical potential in Fig. |3|. In effect here 
is a separation between the phenomenology of quarks of different colours. Quarks 
orthogonal in colour to the diquark condensate (colours 1 and 2) do not conspire to 
break chiral symmetry (so g\ = 0), and rather condense in colour-asymmetric pairs. 
This leaves the quarks colour-parallel (colour 3) to the diquark 3 to spontaneously 
break chiral symmetry in the usual manner (g 2 7^ 0). Thus they acquire an effetive 
mass due to self-interaction and bestow a mass upon the other two quarks through 
a similiar interaction. Since the SU(3) colour symmetry has been broken, these 
masses need not be equivalent, and in fact one finds that 

M 1 = -g 2 , M 2 = 2g 2 . (31) 

In the vacuum, M 1 = 220 MeV and M 2 = 260 MeV. As the chemical potential 
increases, the coexisting diquark and chiral condensates compete for the coupling 
strength; this is the condition of eq. fl2H|) . We find that at relatively low /i, around 
80 MeV, this limited resource becomes entirely consumed by the diquark. This 
is demonstrated in Fig. |3], where one notes the chiral condensate in the colour-3 
channel abruptly vanishes and the diquark condensate assumes its value of phase 
(2). Thus this solution is not only strictly a local minimum, but is also short-lived 
in density. 

The remaining thermodynamic competition is thus between the two phase of 
pure chiral or colour symmetry breaking. For three colours, the critical ratio of the 
condensates (|30| ) is \^3. At /i = we find f/g= 165 MeV/65 MeV > v^- It means 
that at low values of \x the coupling constant A is smaller in the chiral broken phase, 
hence this phase is energetically preferred. Furthermore, this is the A which leads to 
a reasonable effective quark mass of M = 346 MeV, as well as a chiral condensate 
of (ipip) = —(255 MeV) 3 . We thus recover the standard chiral symmetry breaking 
at low baryon density. 

With increasing chemical potential both / and g condensates in their respective 
phases are reduced, however as /i surpasses the effective quark mass M around 300 
MeV, the chiral-breaking g stabilizes while the colour-breaking / continues to de- 
crease in phase (2). At \l ~ 340 MeV condition (|30D is reached, and therefore for 
this and greater chemical potential chiral symmetry is restored and the 'colour su- 
perconducting' phase is realized. At the transition point, the quark gap is A = 115 
MeV and slowly decreasing with rising density. Numerical results for the charac- 
teristic quantities of each phase are plotted in Fig. |]. It is noteworthy that the 
chiral condensate, (ipip) , is distinct from the effective quark mass, M, in that while 
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Figure 4: Condensates for N c = 3 as a function of [i. Shown are the effective quark 
mass M, the quark condensate — (■?/>■?/>) 1,/3 , and the diquark energy gap per quark A. 

the latter decreases with \x the former remains practically at its vacuum value. For 
ji> n c they both vanish. 

7.2 Properties of Quark Matter 

The chemical potential is a useful concept only in that it constrains the particle 
number or density, the physical quantity we actually need to compute. The four- 
fermion interactions we are dealing with are somewhat unusual in that they are 
dependent on the energy a consequence of using the zero mode Green functions as 
the seed of the interaction. 

For this reason we carefully define the density as 

n = J d 4 x ji{x) , (32) 

and note that in this case one does not have the simple form = ip^^ip for the 
current which carries conserved particle number. In our notation this corresponds 
to the current of the high-momentum quarks, which without the zero modes would 
imply that the correct normalization is not preserved. We find the applicable current 
is 
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d 



2/32P2 



+ (L «-> i?) 



(33) 



where the bracketed momentum derivative is understood to be a sum of four terms, 
each with one formfactor differentiated with respect to its momentum argument. 

Although there has been some rearrangement, this is essentially the vertex ([Tj]) 
in which the formfactors have been differentiated with respect to the momentum q^. 
One can verify this current is conserved by considering the Dirac equation, written 
here for the left-handed, flavour- 1 component: 



A 



d A p 2 d^k\d A k 2 



(2< 



12 



{2n)%q + p 2 -k 1 -k 2 



»2 £ k 2 l2TTa2TT}02 
k 2 t U h U $2 e «202- 



^(k 2 ,-^ 2 L ^(k 2 )=0. (34) 



Acting on this with q fJ ,4'Li( y q)-^-, and with similar operations on similar equations 
for the other components, one recovers the required q^j^q) = 0. 

For three colours we have found two phases which specify the stable equilibria, 
each having one broken symmetry. Phase (1) denotes spontaneously broken chiral 
symmetry, and it competes with phase (2) which breaks colour symmetry. As an 
intermediate step in calculating the density, we first examine the quark occupation 
numbers for the two states, n(p) where in this context p = \p\, which is obtained 
via a pi integration and precisely defined as: 



n 



—i 



d 3 p 



n(p) . 



(35) 



With this definition we have not separated the quark and anti-quark distribution 
functions, since the phases we analyse do not necessarily retain these as basic degrees 
of freedom. 

The conserved current written above is a sum of two parts. The first is computed 
with the Green functions obtained through self-consistent solution of the Schwinger- 
Dyson-Gorkov equations, and to the density it contributes 

(-4>KphMp)) = -"' f dp ™- ro/ ^" 1 

r P4 + i[i 



—i 



(2vr) 4 
d 4 p 



(27T) 



(p + i^y 



[2Z l (p) + (N c -2)Z 2 (p)}(36) 



20 



The scalar functions Zi^ip) are complex, and for a given solution set (M\,M 2 , A) 
become 



Zxip) 
Z 2 (p) 



l+a*(p,fi)M 2 -4(3(p,n)A 2 



[1 + a(p, n)M?} [1 + a*(p, /i)M|] - 16/% /i) 2 A 4 
1 



l + a(p,n)Mi ' 



(37) 



The functions a(p, /i) and (3(p,fi), as defined previously, arise as combinations of 
the formfactors evaluated with appropriate momentum arguments. In general, the 
first is complex while the second is purely real and always positive, a consequence 
of this being the expression which accompanies the diquark propagation. Through 
pairing quarks, the complications due to the non-hermiticity of the single-quark 
Dirac operator are avoided. 

In the current's second term the formfactors have been differentiated with respect 
to p 4 . The expectation values were computed to the same O^X 1 ) as followed for the 
Green functions, and this generates another sibling ^(p, /j,) in our family of functions. 



lip, fi) = (j> + i^) a {p + ip,) a (pp(p, n)(pp(j>, fJ>)<P\ 
Combining the two pieces, the density reads 



dp A 



(3* 



n 



-iN f 



d A p 



Pa + ijJ, 



(p + in) 2 
Pa + ifi 



[2Z 1 {p) + {N c -2)Z 2 {p)\ 



■a(p) +7(p) 



^- W )G 1 {p)9i 



+ {N C - 2) 2 G 2 (p)g 2 + (2N C - 5 + J- ) (O ,«/»//, i G,( P ) Ui 



(39) 



The functions Gi^{p) can be deduced from Eqs. (]17D and (p7|). The integrand 
here is complex, however under p 4 reflection we note that a(— Pa) = u*{pa) and 
7(— Pa) = — J*{pa)- This follows from our preserving the non-hermitian contributions 
exactly, again we note that upon p 4 integrate we recover purely real occupation 
numbers and hence quark density. 

We now consider the two phases individually. When phase (1) is assumed, and 
A = 0, the density expression (|39|) simplifies to a form where the integrand contains 
a denomenator of the form (p + ifi) 2 + M(p,fi) 2 . A chiral effective mass is clearly 
present, which upon evaluation of eqs. (pl],p4[) satisfies the same gap equation as 



in the Nf = 1 case, eq. flT3|) f|. The effective mass shifts the Fermi surface, which 



^In fact the same gap equation applies to order A 1 for any number of flavours |27 
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H = 318 MeV 
\i = 360 MeV 
H = 600 MeV 
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Figure 5: Occupation number n(p) vs. p in phase (1), where /i = 318, 360, and 600 
MeV. The corresponding p F = 45, 170, and 450 MeV, respectively. 



for free quarks is at pf = fi, to the solution of pf = \J — A4(pf, aO 2 - This, as 
shown in Fig. |5|, is the only modification of the standard fermionic step function (of 
magnitude 2NfN c = 12) due to chiral symmetry breaking. 

The occupation numbers are more interesting when the quarks fall into the 
colour superconducting state. As is known in BCS theory, Cooper pairing of quarks 
spreads the Fermi surface as fermions pair with momenta above and below the 
chemical potential. As demonstrated in Fig. || this is precisely the result for the 
two quark colours which participate in the diquark. Quarks of the parallel colour 
are unaffected by the energy gap, and furthermore since chiral symmetry has been 
completely restored appear as massless fermions with a degeneracy of 2Np- This is 
the origin of the discontinuity at p = //, since the Fermi momentum remains at this 
point for these free quarks. Thus the plot should be considered a combination of a 
smooth distribution of Bose-condensed diquarks added to a simple step function of 
magnitude 4. To be explicit, eq. (|39| ) can be rewritten as 



n 



(2) 



2N f (N c 



6tt 2 
+2N F 



2) 3 



KpI + v 2 - p 2 ) 



d 4 p 

(2^R+(p-^) 2 M + (p + ^ 



2^21 



l + 4/3(p)A 2 



(40) 



The first term is the density for N c — 2 free quarks, the second a contribution from 
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Figure 6: Occupation number n(p) vs. p in phase (2), where /i = 60, 360, and 600 
MeV. The plot on the left is for all three colours together, the right for the two 
which participate in the bosonic diquark condensate. 

the remaining degrees of freedom which are no longer properly described as single 
quarks. 

Performing the integration over all momenta, we arrive at the density as a func- 
tion of chemical potential. In Fig. [j] this is shown for the equilibria states, demon- 
strating the large discontinuity at the phase transition, where the dashed lines are 
the metastable continuations in each phase. Since the constituent mass of a quark 
in the chiral breaking phase remains near 300 MeV, no density is amassed below 
this chemical potential. This is consistent with our finding a virtually constant 
chiral condensate below the critical chemical potential, in that the absence of phys- 
ical quark density should leave the physical condensate unchanged from its vacuum 
value. Soon after the point at which the Fermi surface moves from zero and physical 
states begin to be filled, the onset of the phase transition has been reached. Thus, 
for pure quark matter in this simple mean- field treatment, the maximum density 
of a purely chiral-broken phase is a rather small n = 0.062 fm~ 3 . Taking a naive 
saturation density for quark matter as three times that for nuclear matter, we have 
n = 0.445 fm -3 , which lies within the coexistence line between the two phases. It 
is only above a density of n = 1.05 fm -3 = 2.53n that one finds the pure supercon- 
ducing phase. For reference, the density of free quarks at this fi would be 2.25no, 
somewhat lower than that which is found in the pairing scenario where the states 
are in some sense compressed through bosonic condensation. 

This extremely low density for the onset of chiral restoration suggests that in the 
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Figure 7: Quark density as a function of chemical potential. The solid line de- 
notes the the equilibrium state of quark matter, the dashed lines the metastable 
continuations. Naive quark matter saturation density lies at ~ 0.45 fm -3 . 

core of heavy nuclei quarks exist in some transitional, 'boiling' state . However, 
we have disregarded certain effects which could modify these results. Considering 
an extended mean-field ansatz of more than the three channels we consider would 
likely result in additional interquark forces, which could lead to enhanced repulsion 
and extend the stability of the chiral-breaking phase to higher density. We have also 
not taken into account the binding energy of the quarks within a nucleon, which 
would not modify our results at the quark level but could lead to some differences 
between quark and nuclear matter. 

Another possibility is to vary the parameters of our theory. In particular, one 
can change the instanton packing fraction which has been fixed at p/R = 1/3 in 
all calculations described here. This in turn scales the strength of the induced 
interaction. However, taking values from half to twice our choice does not lead 
to qualitative differences in the density results; in every case chiral symmetry is 
restored at densities which are a small fraction of that for saturation. 

Finally, we note that the strong prediction one can make from this and similar 
treatments is the presence of a first-order transition which restores chiral symme- 
try. It is a conclusion which is shared by other recent works which consider chiral 
symmetry breaking at the quark level. In particular we note that the density dis- 
continuity which we find at this transition, An m 2.4no, is close to that found in 
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the NJL-type treatment of Berges and Rajagopal [18j (An « 1.5no — 2no) as well 
as that from the chiral random matrix model (An rs 2.5no) ||29|| . 

8 Conclusions 

We have formulated the effective low-energy action for two light fermions induced by 
instantons at nonzero chemical potential. In the resulting fermion vertex we have 
retained the full dependence on both momentum and chemical potential, which 
arises from the would-be zero modes. The overall interaction strength is given by 
a dynamically-determined coupling, which is a consequence of the fact that the 
instanton weight itself is proportional to the fermion determinant. In these respects 
we differ from other studies such as the random matrix model, Nambu-Jona-Lasinio 
models, and alternative instanton approaches. 

In general, introduction of the fermion chemical potential leads to complications 
due to the resulting non-hermiticity of the Dirac operator. This is particularly 
pronounced through the complex nature of the quark interaction, induced by the 
zero-mode solution of the complex eigenvalue problem. By retaining the exact func- 
tional forms without simplifications we are able to avoid any complications from 
individual imaginary eigenvalues. This is an advantage over various numerical tech- 
niques, and should persist for all orders in the instanton density. The formalism 
outlined in this paper allows one to make such an expansion in a systematic way. 

The effective action leads to a competition between two phases, one of chiral 
symmetry breaking and another characterized by diquark condensation. It was 
studied by solving a coupled system of gap equations to first order in the instanton 
density. For two massless flavours and three colours, a simple ansatz allowed us 
a detailed study of the various types of symmetry breaking. Although we did re- 
cover an interesting case of mixed condensation which broke both chiral and colour 
symmetries, it was found to be thermodynamically disfavoured. Considering the 
remaining possibilities, we find that spontaneously broken chiral symmetry is re- 
stored through a first order phase transition, replaced by colour breaking due to 
the formation of a diquark condensate. With our 'standard' choice of the instanton 
ensemble, N/V = lfm -4 and p/R— 1/3, we find the critical chemical potential is 
~ 340 MeV, where the superconducting energy gap is 115 MeV and decreasing with 
rising /x. This translates into the onset of chiral restoration at very low density of 
quark matter. Indeed, the critical \i c = 340 MeV corresponds to the quark density 
of 0.14 no in the chiral broken phase and to the density of 2.53 no in the super- 
conducting phase, where no = 0.445 fm 3 is the quark density corresponding to the 
standard nuclear matter. Taken literally, it suggests that the deep interior of heavy 
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nuclei are in a 'boiling' mixture of the two phases, as suggested earlier in ref. [[16 
It should be kept in mind however, that our calculations were to only first order in 
the instanton density and that we neglected the binding of constituent quarks into 
nucleons. Both circumstances might shift the critical density somewhat. 

Restricting our discussion to the quark level, we computed the fermion occupa- 
tion numbers in the two distinct phases which explicitly demonstrate the physics in 
each system. In the chiral broken phase, the reduced Fermi radius illustrates the 
effective mass, whereas in the diquark phase BCS-type behaviour is demonstrated 
by the lack of a sharp Fermi surface. 
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APPENDIX. Fourier Transforms of Fermion Zero 
Modes 

The use of the exact fermion zero modes in the momentum space tremendously 
simplifies all calculations. The starting point is the exact fermion zero mode in the 
field of one (anti) instanton in x space [ [2~T| which we cite for arbitrary instanton 
position z, size p and orientation given by rectangular N c x 2 matrix U. In the 
spinor representation for the Dirac matrices the zero modes are two-component 
Weyl spinors which can be written as: 



Il(x - z) 

The single instanton solution is apparent in the functions 

P 2 



1 !(.'•) 1 I -, 
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The Fourier transform is defined as 
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and has the structure 

Lorentz symmetry is broken at finite //, and the components become 
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+(p + yu + ip 4 ) [(2p4 + - «(P + A* + ^O/W] } , 
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where the scalar p = |pj, the spatial i = 1 ... 3, and the functions 

Ji(,z±)ifo(^)-/o(^)ifi(^) /i(*±)tfi(;z±) 

Jl± — 5 J2± — — 



Z± 



Z ± 



(46) 



are evaluated at z± = \p\jp\ + (p ± yu) 2 . With these expressions it is explicitly 
verified that the normalization condition holds for any [i: 

" '" $ r (p, Ai) = / [<pl(-fi)<pM + ■ <p(j*)] ■ (47) 
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